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this dispersive effect with the electronic Kerr effect in photonic crystals, and 

attribute the the absorptive features to nondegenerate two photon absorption. 
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1. Introduction 

The interest to optically switch photonic structures has recently gathered momentum due 
to the inherent fastness of the process. While the switching speed of conventional transistors 
is ultimately limited by heat dissipation, no such limitation exists for optical systems in 
the absence of absorption. To date, the dispersion of photo-induced free carriers is widely 
exploited to induce ultrafast changes in the photonic structure's optical properties. Then, 
propagation switching, optical modulation, trapping and releasing photons, frequency- and 
bandwidth conversion, and even ultrafast switching of the density of states are possible [1-6]. 

In the past, several groups have switched Si photonic crystals via free carrier excitation 
[7-9]. This switching mechanism is popular because of its large attainable change in refractive 
index and its high possible repetition rate. Indeed, currently many efforts are devoted to 
decreasing the timescales of switching in polycrystalline silicon (p-Si), achievable through 
implanting additional recombination centers such as ions [10-12]. However, the possible 
recombination time is still limited to several ps. 

Electronic Kerr switching offers the potential to switch instantaneously, with a repetition 
rate limited only by the duration of the pump pulse. The duration of the pump and probe 
pulses then controls the repetition rates, instead of the recombination times of the free 
carriers. Kerr switching could potentially increase the repetition rate from GHz to THz. 
Unfortunately, the intensities required are typically much higher compared to free carrier 
switching [13-15]. Recently, several groups have claimed to have experimental evidence of this 
Kerr nonlinearity in photonic crystals [16-20]. The reasoning in all cases is the instantaneous 
nature of the effect, resulting in a decrease in reflectivity or transmission at one probe 
frequency which coincides with the cross-correlation of the pump and probe pulses. Here, we 
critically evaluate these optical processes at instantaneous (fs) timescales. At pump-probe 
coincidence, we will show that two competing processes change the optical properties as 
witnessed from reflectivity measurements performed over a large bandwidth. From these 
two processes, we propose a general nondegenerate figure of merit (NFOM), which has wide 



implications for all-optical fs switching, ultimately extending to DOS switching and the 
dynamic control of spontaneous emission. 

2. Experimental setup and sample 

Our setup consists of a regeneratively amplified pulsed Ti:Sapph laser (Spectra Physics 
Hurricane) which drives two independently tunable optical parametric amplifiers (OPAs, 
Topas). The OPAs have a continuously tunable output frequency between 0.44 and 2.4 eV, 
with nearly bandwidth limited pulses (relative width of 1.33%, bandwidth limited equivalent 
to rp = 110 fs) with pulse durations of rp = 140 ± 10 fs (measured at -Ep ump = 0.95 eV) and 
a pulse energy of at least 20 /iJ. The probe beam is normally incident 9 = 0° on the sample, 
and is focused to a Gaussian spot of 32 /im FWHM (at -Enrobe = 1-24 eV) at a small angular 
divergence NA = 0.02. The probe intensity Ip ro be of 3 ± 2 GWcm" 2 was around 10 times 
lower than the pump intensity Ip ump to ensure that nonlinear effects caused by the probe 
were negligible. The reflectivity was calibrated by referencing to a gold mirror. A versatile 
measurement scheme was developed to subtract the pump background from the probe signal, 
and to compensate for possible pulse-to-pulse variations in the output of our laser [21,22]. 

Figure [1] shows a high resolution scanning electron micrograph of the photonic woodpile 
structure made by Jim Fleming at Sandia National Laboratories [23]. It consists of 5 layers 
of p-Si rods (n' = 3.45) stacked orthogonally upon one another, the n-th layer shifted by half 
an interrod distance with respect to the (n+2)th layer, with a total thickness of L = 780 nm. 
This structure gives rise to a diamond lattice, for which a first-order band gap is predicted 
[24] . The dimensions of the rods were chosen so as to aim the center of the band gap around 
telecom frequencies (-Ereie = 0.735 eV). The electric field of the probe beam is polarized 
along the [110] direction of the crystal, that is, perpendicular to the first row of rods. The 
last row of rods is supported by a 70 nm thick SiN layer of refractive index 2. 

The sample has been characterized extensively elsewhere [22]. Over the sample, the mutual 
alignment of the rods differs. These differences are spatially separated, and allow the sample 



to be divided into 16 domains, which we denote by Al through D4 in analogy to a chessboard. 
Here, we perform measurements on both Al and D4, and are thus able to validate our results 
for different sample conditions. The symmetry of D4 is face centered orthorhombic, while 
Al is body centered orthorhombic. 

3. Linear reflectivity 

Figure [2] shows a linear reflectivity spectrum of the woodpile photonic crystal. The high 
peak centered at Eq = 0.9 eV corresponds to the T — X stop gap in the band structure, 
which is part of the 3D photonic band gap of Si woodpile photonic crystals [21,25,26]. The 
maximum reflectivity of 95 ± 2% and the broad width of AE = 0.46 eV FWHM (full width 
at half maximum) confirm the strong interaction of the crystal with light. The interaction 
strength S is characterized by the relative gap width S = AE/E [27]. We note that the 
large photonic strength of S = 47% is the largest measured for any photonic crystal so far. 
We can associate a length scale with the photonic strength. The Bragg length Lb is the 
length at which the incoming intensity has dropped to 1/e of its initial value, and is given 
by [28] 

where A is the wavelength in the structure of the stopgap. We deduce a Bragg length of 
Lb = 840 nm, corresponding to 1.06 unit cells. In contrast, other 3D structures with a lower 
photonic strength have a much higher Lb'- For example, TiC>2 inverse opals have Lb ~ 3.3 /mi, 
polystyrene opals have Lb ~ 5.1 /iin and SiC>2 opals Lb ~ 6 /tm [28]. 

The high quality of the sample also causes the steep edges of the peak, which we will in first 
instance probe at the blue edge of the stopband. At these frequencies, all induced changes 
in optical properties of the Si backbone will cause a large change in reflectivity, important 
for switching applications. The probe frequencies (indicated by the shaded box) span both 
the blue edge of the stopband as well as the fringes. These have been tentatively assigned 
to Fabry- Perot-type interferences [29], but the interpretation of the reflectivity spectrum of 



woodpiles is still an open question [22]. 

4. Switched reflectivity vs. time at one probe frequency 

The pump frequencies Ep ump were chosen as to tune through half the electronic band gap of 
silicon, where the gap is Eq = 1.12 eV. Both the Kerr coefficient n<i and the degenerate two- 
photon absorption coefficient (3u have recently been shown to vary strongly in the vicinity of 
\Eq [13-15,30-32]. On domain Al we measured differential reflectivity AR/R vs. probe delay 
At as a function of Ep ump (figure [3]). On this domain, the probe frequency -EWobc = 1-13 eV 
corresponds to the foot of the blue stopband edge. At positive probe delays At = 1 ps and at 
pump frequencies -Ep ump > \Eq = 0.56 eV, we observe a large positive differential reflectivity 
(AR/ R)fq due to a blueshift of the photonic features. This shift is caused by a decrease of 
the refractive index of the p-Si backbone that agrees well with a Drude description of the 
excited free carriers [21]. When Ep ump is reduced to below \Eq, the dispersion vanishes, 
since no free carriers are excited by a two-photon process. 

At pump-probe coincidence, near At = ps, we observe a trough of magnitude 
(AR/R) coinc . The width of this trough varies between 240 fs (-Ep ump = 0.52 eV) and 160 fs 
(-Ep ump = 0.62 eV). These widths agree reasonably well with the expected cross-correlation 
duration for two pulses (rp x \[2 = 200 fs). As Ep ump increases to above ^Eq, the magnitude 
of the instantaneous feature (Ai?/i?) coinc grows until it saturates at -Ep U m P = 0.62 eV, just 
over \Eq. In all cases, instantaneous switching features are observed. 

5. Switched reflectivity vs. frequency 

5. A. Reflectivity vs. frequency at coincidence 

In order to investigate the precise nature of (Ai?/i?) coinc , we have measured the differential 
reflectivity at coincidence on domain D4 for a wide range of probe frequencies, which we 
plot in figure 11(c). We make three important observations. Firstly, with increasing £p ump , 
the variation in the data becomes stronger. Secondly, the differential reflectivity becomes 



increasingly negative with both increasing pump and increasing probe frequency. These neg- 
ative values indicate optical absorption. Thirdly, at 1.22 eV, a peak appears in the differential 
reflectivity that corresponds to a red edge in the linear reflectivity (figure IH(a)); likewise, a 
trough appears at 1.18 eV which corresponds to a blue edge. This clear dispersive behavior 
is indicative of a redshift of the photonic features related to a positive change in n! . Thus we 
conclude that when pump and probe are coincident, two effects contribute to the differen- 
tial reflectivity: Nonlinear dispersion, which we attribute to the electronic Kerr effect, and 
absorption. We conclude that the negative differential reflectivity is due to an instantaneous 
nondegenerate two photon process, since this absorption is seen to increase for an increase of 
either pump or probe frequencies, and since the sum of pump and probe frequencies exceeds 
the absorption edge of p-Si. Recent work [16-20] has attributed the behavior at coincidence 
to the Kerr effect exclusively, even though dispersive data are lacking. Here, for the first 
time we identify the different instantaneous contributions to all-optical photonic switching. 
Therefore, this also represents the first identification of Kerr switching in photonic crystals. 

5.B. Reflectivity vs. frequency at At = 1 ps 

Before quantifying the two different instantaneous contributions from the data, we will first 
briefly describe the measured effects of the free carriers on the spectral properties. Figure 
H(b) shows the differential reflectivity at a positive probe delay of At = 1 ps that is caused 
by free carriers. Because the differential reflectivity varies around 0, we conclude that the 
induced absorption is minimal. The peak at 1.2 eV corresponds to the blue edge of the linear 
reflectivity peak at 1.175 eV, the two troughs at 1.23 eV and 1.16 eV corresponds to the 
red edges of the reflectivity peaks at 1.175 and 1.25 eV, respectively. Therefore, the data are 
consistent with a blueshift, or a decrease of the refractive index An' as a result of the free 
carriers. 



5. C. Model: The Extended Scalar Wave Approximation 

To quantify the two nonlinear contributions at coincidence as well as the free carrier effects, 
a model of the woodpile's linear and nonlinear spectral properties is necessary. To predict the 
change in complex refractive index, we require a physical model which is subject to two re- 
quirements: First, this model must reproduce the linear reflectivity Rq over a given frequency 
range. This requirement is necessary because the measured differential reflectivity, defined as 

^ - = ' R ~ ° , depends sensitively on R . Second, the calculated spectral properties must 
scale correctly with the induced complex n. 

Here, we employ a heuristic model to calculate both the linear as well as the dynamic 
reflectivity, which is an extension of the scalar wave approximation (SWA) [33]. The dielectric 
function in beam direction e{z) and the lattice spacing can be written in terms of their 
Fourier components Uq and G, respectively. For binary structures composed of two different 
dielectric materials, it can be shown that [28] 

U G = Aef G ,G^0, (2) 

where Ae is the difference of the dielectric constants, and fo is the Fourier Transform of the 
indicator function f(z), where e(z) = ei + Aef(z). From a change in Ae, a new Uq can readily 
be obtained. The electric field inside the crystal can now be calculated by considering only 
the first two bands, k (incident) and k — G (Bragg diffracted). Outside the crystal, the field 
can be determined from the two boundary conditions the two bands are subject to. Thus, the 
transmission and reflectivity can be determined. Recently, Euser et al. employed an exact 
modal method (EMM) to calculate the broadband reflectivity of the woodpile photonic 
crystal [25]. Although the agreement was good for low frequencies, the reflectivity could not 
be reproduced faithfully beyond the stopgap [21]. To obtain a better agreement, the number 
of plane waves would have had to have been increased, for which the calculation diverges 
slower. Also, Deubel et al. acheived good agreement at the blue edge with a scattering matrix 



model, although some salient features were still lacking [34]. We verified that at the relevant 
probe frequencies, the induced shift obtained with SWA matches the EMM (see appendix 
[A}. 

We fit the linear reflectivity obtained by the scalar wave approximation to the region given 
in fig. S(a). The best agreement to the measured reflectivity is obtained with G = 27r/(393 
nm) and Uq = —0.8. The calculated spectral feature shown corresponds to the first two 
Fabry-Perot fringes at the blue edge of the stopgap. In the model, we included the dispersion 
and absorption of p-Si (from [35], see fig. [5]), and have taken bulk Si as the substrate. We 
verified that the thin layer of SiN (n' — 2, d — 70nm) had no effect on our analysis (appendix 
IB"]) . Since the optical features in the reflectivity spectrum are much broader than the probe 
bandwidth, no convolution of the calculated spectrum with the probe is necessary. 

5.D. Interpretation of spectra at coincidence 

We match the calculated dispersive differential reflectivity to the strong signature of a 
redshift at -Ep ro bc = 1.2 eV by fitting two independent parameters An' = 2ra 2 /pump and 
n" = Ao/3i2-fpump/27r over the range 1.18 eV < -Eprobc < 1-23 eV by minimizing the least 
squares. Here, j3u is the nondegenerate two-photon absorption coefficient and Ao the free 
space wavelength. The differential reflectivity was obtained as follows: using the G obtained 
previously from the fit to the linear spectrum, the Ug was changed via the changed Ae 
(equation [2]), and the switched reflectivity spectrum was calculated from the new Uq. Error 
margins for the parameters were taken as the values at which the least squares values had 
changed by 10%. 

Figure H^c) shows the fits to two spectra obtained at two different -Ep ump . The remarkable 
agreement of this approximate model with the data confirms that the subtle interplay of 
nonlinear dispersion as well as nondegenerate absorption leads to the observed features. At 
frequencies below 1.2 eV, the fitted curves tend to be lower than the measurements. Above 
1.2 eV, the calculated differential reflectivity is mostly higher than the measurements. These 



deviations are due to our fit not taking into account the frequency dependence of /3 12 , but 
setting it to one value for all probe frequencies. 

Figure [5](a) and (b) shows the absorptive and dispersive part of the induced refractive 
index at coincidence, respectively, as returned from the fits. The value at Ep nmp = 0.59 eV 
is excluded due to an experimental artifact in that measurement. Figure [5](a) shows that 
the imaginary refractive index n" increases for pump frequencies -Ep ump > 0.55 eV. These 
frequencies correspond to a sum of pump and probe frequencies E Pump + E Prohc > 1.75 eV. 
Since these values are larger than the optical gap of p-Si (1.4 eV < E opt < 1.6 eV, depending 
on the degree of polycrystallinity [36]), it is clear that the simultaneous presence of pump 
and probe photons induces optical absorption. The resulting values of the nondegenerate two 
photon absorption coefficient are /3\2 = 0.11 cmGW -1 at -Ep U m P = 0.54 eV and /?i2 = 0.19 
cmGW" 1 at E Pump = 0.62 eV. 

From recent calculations [37], the iVph-photon process in a direct transition can be viewed 
as a single-photon process with an energy gap rescaled to Eg/N p ^, where N p ^ is the number 
of photons. In other words, (-Ep ump + -Ep ro bc)/2 = E' in the present two beam case, where 
E' is a rescaled energy. Therefore, we can compare the spectral shape of our extracted 
nondegenerate absorption coefficient to the linear absorption. We therefore also plot the 
linear absorption of p-Si [35], similar to our woodpile crystals [23]. The behavior of the 
linear absorption agrees reasonably with that of the nondegenerate absorption, supporting 
calculations from [37]. For our intensities, the nondegenerate absorption remains significantly 
below the linear absorption. Pumping or probing at higher frequencies causes the absorption 
to increase, as do higher pump intensities, but not higher probe intensities [38]. 

The instantaneous change of the real part of the refractive index An' is plotted in figure 
Mjo). Again the value at E Pnmp = 0.59 eV is excluded from the discussion. At low pump 
frequencies of 0.52 eV, the change is 0.44 x 10~ 3 , and it more than doubles to 1.1 x 10~ 3 at 
-Epump = 0.56 eV. At Ep ump = 0.62 eV, An' decreases slightly. The resulting values of the 
Kerr coefficients are n 2 = 7.3 x 10~ 6 cn^GW" 1 at E Pump = 0.54 eV and n 2 = 12.7 x 10~ 6 
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cm 2 GW~ 1 at £p ump = 0.62 eV. We compare our measurements to recent measurements 
of the degenerate nonlinear dispersion of bulk Si [13-15]. In order to be able to make the 
comparison, we made use of the rescaling condition. We find excellent agreement to the 
data from [13, 15], but surprisingly the measurements of [14] are an order of magnitude 
higher, despite all being measured with the same z-scan technique. We also compare our 
measured data to a relation for n 2 derived for direct-gap semiconductors [39] multiplied by 
our pump intensity ip ump . The functional form agrees well with our data. At higher probe 
frequencies, the theoretical relation differs from the data of Refs. [13,15]. The agreement can 
be improved by employing a theory for the degenerate nonlinear properties of an indirect 
semiconductor [37] using the rescaling condition [14]. 

Having identified the two contributions at pump-probe coincidence, we briefly discuss 
recent photonic crystal switching experiments in the light of our new interpretation. Becker 
et al. performed time-resolved non-degenerate pump-probe spectroscopy on a-Si and p-Si 
inverse opals [8]. In their data on the a-Si sample, sharp features indicative of instantaneous 
absorption are apparent for E Pump = 0.674 eV and -Ep ro bc > 0.653 eV, or ©Total > 1-327 eV 
at pump-probe coincidence. Their p-Si sample reveals these sharp features at E Pump = 0.717 
eV and E Pmhe > 0.539 eV (©Total > 1-256 eV). Because of the lower optical gap E Qpt for 
p-Si with respect to a-Si, we expect nondegenerate absorption to take place at a lower ©Total 
than for the a-Si sample, in agreement with their data. As a second example, we discuss 
the contribution by Hache and Bourgeois [16]. Here, the authors probe a Si/SiC>2 Bragg 
stack at the red edge of the stopgap. At coincidence, they observe a significant increase in 
transmission. The observed feature can only be explained by the breakdown of destructive 
interference due to nondegenerate absorption for ©p ump = 0.725 eV and E Prohc = 0.821 
eV, giving .©Total = 1-546 eV, comparable to Eo pt of a-Si. The Kerr effect without induced 
absorption would lead to a red shift of the stopgap and thus to a decrease in transmission. 
In a similar manner, the contributions by Mazurenko et al., Hastings et al, and Mondia et 
al. misinterpret the sharp trough in differential reflectivity as the Kerr effect [17-19]. 
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5.E. Interpretation of ps switching 

Using the same extended SWA from section 15.01 we can also interpret the picosecond free 
carrier behavior. Therefore, the calculated differential reflectivity can be fitted to (AR/ R)-pc 
by calculating the complex refractive index change with the Drude model (figure H^b)). The 
carrier density N is given by 

2-Ep ump e 
where fin is the degenerate two photon absorption coefficient. At our pump frequencies, 
we can neglect linear absorption. For carrier densities N < 10 28 m~ 3 , well above densities 
created in our experiments, the change in dielectric constant is [40] 

A6 fc = -(^) 2 (l-»(«7i 3 )- 1 ) ) (4) 

where uj 2 p = Ne 2 /(m*eo) is the plasma frequency squared in (rad/s) 2 and e and m* are the 
charge and effective optical mass of the carriers, and uj the probe frequency in rad/s. Equation 
2] is valid for (cut^) -1 « 1, which is the case for our probe frequencies and damping times. 
For a given 7p ump therefore, we can fit the differential reflectivity resulting from a changing 
Pn to (AR/R) FC - 

The /?n versus E Purap as returned from the fits is plotted in figure [BJ For E Pump < \Eq, 
j3n tends to 0, as expected. Increasing -Ep ump to above 0.54 eV, there is a marked increase. 
At -Epump = 0.62 eV, (3\\ has increased to as much as 5 ± 1 cmGW -1 . This value compares 
favorably to other values of p-Si samples: values as high as j3n = 60 ± 15 cmGW -1 have been 
reported [22]. Indeed, a systematic study by Ikeda et al. showed the nonlinear coefficient 
increases with increasing amorphousness of the sample, with (3n = 4 cmGW" 1 for p-Si 
samples, and up to 120 cmGW -1 for amorphous silicon [41]. To summarize this comparison, 
we plot our extracted values of /?n along with j3 u of both p-Si and c-Si sample in fig. [6j 

We are now able to calculate crucial parameters to characterize the quality of the switching 
process in terms of absorption and homogeneity. We find that the carrier induced absorption 
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length is £ a b s = 123 ± 50 yU.ni, much longer than the Bragg length Lb, confirming that the 
absorption is negligible indeed. The pump homogeneity length that describes the spatial 
homogeneity of carriers is £hom = (-fpump/Sii)^ 1 = 44 ± 10 /iin, or 170 ± 40 Lb at Ep nmp = 
0.62 eV and 7p ump = 46 ± 5 GWcm" 2 . This calculation highlights the advantage of two- 
photon absorption [30] to one-photon absorption in which absorption strongly limits the 
pump homogeneity [17,18,42]. We conclude that for our switching conditions, both the pump 
and the probe absorption length remain well above L, and that thus instantaneous switching 
yields favorable conditions for our pump and probe frequencies, and pump intensities. 

6. Interpretation of switched reflectivity at one probe frequency 

Having analyzed spectra of the reflectivity at coincidence and at At = 1 ps, we return to 
figure [3j Experimental conditions had prevented us from taking more than 4 probe frequency- 
and time-resolved spectra, mostly related to the long measurement times. Having identified 
two competing processes at coincidence, we can now focus on time-resolved spectra which 
depend on pump frequency only. Figure [7| shows (Ai?/i?) coinc as a function of £p ump . Here, 
Ip ump has been corrected for three important pump beam parameters: i., the increase in 
intensity with increasing -Ep ump because of the decreasing Gaussian focus [43]; ii. the change 
in pump intensity due to the frequency dependent reflectivity of the pump, see figure EJ and 
iii. the measured 10% change in pump power over the frequency range of the OPA. For L, 
(AR/R) coinc is taken to vary quadratically with -Ep ump , assuming that (Ai?/i?) coinc changes 
linearly with Jp ump . For ii, (Ai?/i?) coinc was linearly corrected by a factor (1 — _R(£p ump )), 
where -Rp U mp(-Epum P ) is the reflectivity of the pump at the pump frequency. Finally, in iii., we 
corrected (Ai?/i?) coinc linearly with /p um p- The two largest corrections are by 80% (-Epump = 
0.52 eV) and 34% (-Ep ump = 0.56 eV), but otherwise the corrections are less than 20%. 

In figure [7] we observe an increasing negative differential reflectivity — (AR/R) coinc with 
increasing -Ep ump up to 0.62 eV, then a decrease until 0.69 eV before a subsequent increase. 
With the analysis of the previous data (fig. [5]), two important regimes can now be identified: 
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for -Ep um p < 0.69 eV (E Total < 1.82 eV), (AR/R) coinc is mostly governed by dispersion. For 
-Epump > 0.69 eV, (AR/R) coinc increases again due to a combination of both dispersive (n 2 ) 
as well as absorptive parts (flu)- 

In figure [7J we also show the differential reflectivity expected from the extended SWA. 
Because this different reflectivity spectrum was obtained on domain Al, we first have to 
apply the SWA to the linear reflectivity from Al (not shown). We obtain best agreement for 
G = 2?t/(363 nm) and XJq = —1.2. From the theoretical relation for a degenerate 7i2 of a direct 
bandgap semiconductor [39], (A_R/i?) coinc is readily obtained. Since a nondegenerate theory 
for indirect bandgap semiconductors is not yet available, the theory has to be interpreted 
cautiously with respect to the measurement. For Ep nmp < 0.7 eV, the shape of the theoretical 
relation agrees excellently to our measurements. A strong divergence from the data is found 
for -Ep um p > 0.7 eV, or -©Total > 1-83 eV. We interpret this threshold frequency as marking 
a transition from dispersive behavior of (Ai?/i?) coinc to absorptive behavior. 

To verify the consistency of the time- and frequency resolved (AR/R) CO i nc taken on domain 
D4 (see section I5.AJ) with the pump frequency resolved (AR/R) coinc taken on domain Al, 
both ?7,2 and /3u (figure [5]) were inserted into the SWA at Al, and plotted in figure [71 The 
obtained differential reflectivity is lower by a factor 2. Reasonable agreement to the measured 
differential reflectivity is found, from which we conclude that the extended SWA provides 
consistent description of our various data, measured at different pump intensities and on 
different sample domains. 

7. Nondegenerate figure of merit for fs switching 

To quantify how useable a material is for optical switching, it is vital to consider a figure of 
merit (FOM). A high FOM then is tantamount to a large phase shift, desirable for switching 
applications. Previously, Garmire proposed a degenerate FOM, a measure which indicates 
the maximum achievable phase shift of 2iix FOM, given that the device length is equal to 
the pump homogeneity length [44]. The absolute length of the device then depends on the 
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intensity used. From this requirement it follows that, 

FOM(© Pump ) = - n2(Ep 7 } (5) 

-APumpPll {-t^Pump) 

where Ap ump is the pump wavelength. Here we extend this notion to the general case by 
including nondegenerate absorption: Not only is the pump absorbed because it excites free 
carriers, but the probe is absorbed in the presence of the pump. The nondegenerate FOM 
(NFOM) is then [45] 

NFOMfE E ) = n 2(Fp ump ) ,g. 

Apump/3ll(Ep ump ) + Ap ro bc/9l2(Epump, Ep ro bc) 

Because of the much lower probe intensity, we can neglect ^(-Eprobc) in the numerator and 
/ 5 11 (£ , p ro b c ) in the denominator. The NFOM has been calculated from our measured data 
as follows: n 2 (Ep urap ) has been derived from figure 0(b), where we have used the functional 
form of ^(-Epump) from [39], scaled by a constant factor to match the magnitude of our 
data. An analytic form of f3\2 was derived by fitting an exponential to the linear absorption 
coefficient a(©p ump + E Prohc ). The exponential behavior of the F absorption coefficient a 
with -Ep ro bc close to the gap region is well documented [46]. From a(E Pump + E PTobe ), we 
deduce p\ 2 = a(E Pump + E Prohc ) / (2I Pump ) from equation Finally, an analytic expression 
for Pn(Ep ump ) is obtained from the theoretical relation for Pn(E Purap ) for direct bandgap 
semiconductors from [39], scaled by a factor to match our measurements (fig. [6]). 

Figure [S] shows the NFOM for different ©p ump and E Prohe . At total frequencies ©Total = 
-Epump + ©Probe below the optical gap of p-Si (-©Total < ©Opt ~ 1-5 eV), the nondegenerate 
two-photon absorption is small, tantamount to a long probe absorption length, which is 
desired. If in addition ©p um p < \Eg, the pump absorption length becomes infinite. Higher 
order photon absorption has to be considered to quantitatively estimate the NFOM in this 
region of parameter space [47]. Around ©p U m P = 0.8 eV, n 2 changes sign and thus the 
NFOM approaches 0. For all other frequencies, NFOM depends strongly on both pump 
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and probe frequencies. To achieve a high NFOM, it is best for -©Total to remain below the 
optical gap, while -Ep ump should remain just below \Eq- In that case, n-i is high, while both 
fin and f3\2 are kept low. We have indicated this line of maximum NFOM in fig. [HI These 
experimental conditions were recently realized by us on a GaAs/AlAs microcavity [48]. Our 
present measurements have taken place in the parameter space indicated by the box. In this 
region, we find NFOM between 6 x 10~ 4 for -Enrobe = 1-24 eV and -Ep ump = 0.75 eV, and 
NFOM = 0.05 for -Ep ro be = 1-09 eV and -Ep ro be = 0.58 eV. The only other measurement 
of 77-2 on a Si based photonic structure is from Ref. [16]. We conclude that for a judicious 
choice of pump and probe frequency, the NFOM can be maximized by limiting both pump 
and nondegenerate probe absorption. 

8. Conclusion and recommendations 

Using nondegenerate tunable pump-probe spectroscopy on a Si woodpile photonic band 
gap crystal, we have for the first time identified Kerr switching in photonic crystals. For 
pump+probe frequencies greater than the optical gap of p-Si, we additionally observe ab- 
sorptive features, which we identify with nondegenerate two photon absorption. Using a 
heuristic model, we estimate the maximum change in the real part of the refractive index 
to be An' = 1.1 x 10~ 3 , giving a Kerr coefficient of n^ = 17 x 10~ 6 , which is in good 
agreement to values in literature. The maximum imaginary part is 1.1 x 10~ 2 , corresponding 
to /3i2 = 1.19 cmGW" 1 at E Punip = 0.62 eV. This complex refractive index gives rise to 
reflectivity changes of 1%. From the interplay of these two effects we derive a nondegenerate 
figure of merit for instantaneous switching (NFOM). For a total pump+probe frequency in 
excess of the optical gap of Si, the NFOM is low, while opportune switching is achieved for 
low probe frequencies, and for Ep ump ~ 1/2Eq. 

At probe delays of At = 1 ps, we find that both pump and probe absorption is negligible, in 
agreement with a Drude description of the free carriers. The differential reflectivity spectra is 
purely dispersionlike, giving changes in reflectivity of up to AR = 2.5%. From this value we 
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deduce the maximum two-photon absorption coefficient to be flu = 5 cmGW -1 , in agreement 
to values stated in literature for p-Si. 

In spite of the challenges in analyzing nondegenerate optical properties of a photonic band 
gap crystal, we have measured a clear difference in response of -©Total < ©Opt and ©Total > 
©Opt- This difference has ramifications on future experiments on nondegenerate instantaneous 
switching photonic bandgap crystals. DOS switching, and ultimately spontaneous emission 
switching is limited by the magnitude of pump and emission frequencies. 

To exploit Kerr switching, that is to maximize NFOM, following improvements are pro- 
posed. In first instance it seems useful to pump and probe with the same polarization. For 
collinear polarization, a factor 3 improvement is predicted with respect to orthogonal po- 
larization [49] . Furthermore, the optical properties of the photonic crystal could be tailored 
so as to enhance the field. Calculations by us on a ID Distributed Bragg Reflector have 
revealed that the field can significantly be enhanced when pumping at either the red or blue 
stopband edge, while the position of the field maximum depends critically on the exact pump 
frequency. 

While a high-bandgap semiconductor will be less prone to nondegenerate absorption, no, 
scales with Eq 3 , thus reducing the NFOM. The absorption can also be diminished by a lower 
pump intensity, which also goes at the expense of An'. We therefore conclude that judicious 
choice of semiconductor material, and pump and probe frequencies are required to usefully 
exploit Kerr switching. 
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List of Figure Captions 

Fig. 1 High resolution scanning electron micrograph of the surface normal to [001] of a Si 
woodpile crystal at domain D4. The width and thickness of each rod is 175 ± lOnm and 
155 ± lOnm, respectively. 

Fig. 2 Linear reflectivity spectrum of the woodpile photonic crystal measured normal to 
the [001] direction at a sample domain D4 shown in figure [TJ The E- field is perpendicular 
to the first row of rods. A stopgap near 0.9 eV gives rise to a high maximum reflectivity of 
95 ± 2% and has a broad relative width of 47%, indicating a high photonic strength. At high 
frequencies > 1.2 eV, the spectral features are attributed to Fabry-Perot-type fringes. The 
pump frequencies (shaded box) were tuned at the red edge through half the electronic band 
gap Eq = 1-12 eV of silicon (vertical dashed line), and the probe frequencies at the blue edge 
of the stopband. The dashed curve is a calculation with the Scalar Wave Approximation in 
the region of interest. 

Fig. 3 Differential reflectivity AR/R versus probe delay At taken at different pump fre- 
quencies -Ep um p and at probe frequency -Ep ro bc = 1-13 eV, sample domain Al. At At = 
ps, the pump and probe are coincident, and the differential reflectivity (Ai?/i?) coinc de- 
creases. At At = 1 ps, the differential reflectivity (AR/R) FC has increased due to the dis- 
persion of the free carriers (FC). The peak pump intensity varies between ip ump = 10 ± 1 
GWcm- 2 (£ Pump = 0.516 eV) and J Pump = 25±2 GWcm- 2 (£ Pump = 0.75 eV), and the probe 
intensity was Ip ro be = 3 ± 2 GWcm -2 . 

Fig. 4 (a): Measured linear reflectivity versus probe frequency (open squares) on domain 
D4 compared to the scalar wave approximation (solid curve), (b) The differential reflectivity 
(AR/ R)fc at delay At = 1 ps caused by the free carriers shows mostly dispersive features, as 
seen from the symmetric variation of AR/R around 0. When increasing the pump frequency 
from 0.54 (solid triangles) to 0.62 eV (open triangles), these dispersive features increase 
in magnitude, (c): At At = ps, dispersive as well as absorptive features are observed in 
(Ai?/i?) coinc , increasingly so when increasing the pump frequency from -Ep ump = 0.54 eV to 
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0.62 eV. For E Pump = 0.54 eV, we find An' = 0.44 x 10~ 3 and n" = 0.46 x 10" 3 , while for 
-Epump = 0.62 eV the fits give An' = 1.1 x 10~ 3 and n" = 1.1 x 10~ 2 . The pump intensity- 
was J Pump = 36 ± 4 GWcnr 2 (£ Pump = 0.54 eV) and J Pump = 46 ± 5 GWcm- 2 (E Pump = 0.62 
eV). 

Fig. 5 Instantaneous changes of the complex refractive index, (a): Imaginary refractive 
index (n", solid squares) vs. pump frequency -Ep ump (lower axis) obtained from fits in fig. |U 
The dashed curve is a fit of n" to an exponential. The right scale shows the corresponding 
absorption coefficient. For comparison, we plot the absorption coefficient for low pressure 
chemical vapor deposited p-Si (solid curve) annealed at 545°C, similar to the backbone of 
our woodpile crystals. Upper abscissa is the sum of the pump and probe frequencies i?Totai- 
(b): Change in real refractive index An' (left scale, solid squares) vs. pump frequency £p U m P 
from our nondegenerate measurements. The theoretical relation predicted by [39] is also 
shown (solid curve, right scale). We plot degenerate measurements for comparison: Open 
squares from [13] and open circles from [15]. Data from [14] exceeds the scale by one order 
of magnitude. 

Fig. 6 Nonlinear coefficient j3n vs. pump frequency (solid squares) obtained from data as 
in fig. |U Dashed curve is a calculation from Ref. [39] for c-Si. Other data are from [13] (c-Si, 
open squares), [15] (c-Si, open circles), [14] (c-Si, triangles) and [41] (p-Si, diamonds). 

Fig. 7 Differential reflectivity at coincidence vs. pump frequency extracted from figure 
[3] (domain Al) and corrected for pump beam parameters (squares, left scale). Circles are 
complex values of n extracted from fig. [5] (domain D4) and reinserted in the extended SWA at 
domain Al. We have calculated the differential reflectivity expected from a nonlinear, purely 
dispersive direct bandgap material (solid curve) [39] by inserting the theoretical relation for 
ri2 into the extended SWA. The 'purely dispersive assumption' deviates from the data for 
pump frequencies above E Punip = 0.69 eV, or E Total = 1.82 eV. 

Fig. 8 Nondegenerate Figure of Merit (NFOM) for instantaneous switching vs. pump and 
probe frequency for polysilicon (see eq. EJ). For Ep ump < \Eq (=0.56 eV), pump absorption 
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is absent as we do not consider more than second order photon absorption. If in addition 
-Epump + -Eprobc < -^Opt = 1-5 eV, the total absorption is leading to large NFOM outside of 
the scale of graph (crossed area). The black curve indicates the parameter space for which 
a maximum NFOM is reached if three photon absorption is neglected. At Ep nmp = 0.8 eV, 
^(-Epump) crosses cm 2 GW _1 , and so the NFOM vanishes. For E Total » E 0pt , the NFOM 
is < 1 x 10 -5 (hatched area). The parameter space in our present measurements are bounded 
by the box (squares). The only other measurement of n^ on a Si based photonic structure 
is [16] (circle). 

Fig. 9 (a) Relative change in frequency AE/E for a given change in the relative average 
refractive index An av /n av . For comparison, the change according to Bragg's law AE/E = 
—An av /n av . (b) Change in reflectivity for a given change in imaginary part of the refractive 
index n" . C 

Fig. 10 Reflectivity calculated with the scalar wave approximation, with and without the 
thin layer of SiN. 
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A. Comparison of the scalar wave approximation to the exact modal method 

To assess whether the changes in spectral properties resulting from a change in n' and n" 
agree with the induced changes in the EMM, we compare the spectral changes for a given 
n' or n". Figure [9](a) shows the relative change in frequency for a given change in average 
refractive index. The shift AE due to An' at the same probe frequency is compared. We 
find that for changes An' < 0.2, which is more than two orders of magnitude larger than our 
expected An' = 1 x 10" 3 , the agreement is even better than 5%. The good agreement for both 
An' and An" proves indeed that the two requirements are sufficient for the description of 
the dynamic reflectivity. Therefore, spectral properties associated with the vectorial nature 
of the probe, such as TE-TM mixing [52] or polarization effects [29], have a negligible effect 
on the change in reflectivity. In figure [9]Jb) we show the expected change in reflectivity 
for a given change in n". To include absorption in the SWA, the reflectivity is reduced by 
exp(— 2L/£ a bs), where £ a b s = Ao/(47rn"0). The factor of 2 in the exponent accounts for the 
beam passing through the sample twice. For -Ep ro bc = 1.19 eV, and for An" < 0.1, which is 
around 10 times larger than expected, the agreement is better than 5%. 

From the good agreement between the exact modal method (EMM) and the extended 
scalar wave approximation, we conclude that the heuristic accurately predicts the change in 
complex refractive index. 

B. Effect of SiN on linear reflectivity spectrum 

In order to verify that the layer of SiN supporting the woodpile has negligible effect on the 
linear reflectivity, a calculation of the SWA is performed with this layer. Figure [TOl shows the 
reflectivity spectra, without and with the supporting SiN layer (70 nm). The spectra shifts 
by 6.2 meV at the most in our region of interest. Therefore we conclude that the inclusion 
of this layer has only a negligible effect on the linear spectra. 

Figures 
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Fig. 1. High resolution scanning electron micrograph of the surface normal to 
[001] of a Si woodpile crystal at domain D4. The width and thickness of each 
rod is 175 ± lOnm and 155 ± lOnm, respectively. SEMD4.eps 
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Fig. 2. Linear reflectivity spectrum of the woodpile photonic crystal measured 
normal to the [001] direction at a sample domain D4 shown in figure HJ The 
E-field is perpendicular to the first row of rods. A stopgap near 0.9 eV gives 
rise to a high maximum reflectivity of 95 ± 2% and has a broad relative width 
of 47%, indicating a high photonic strength. At high frequencies > 1.2 eV, the 
spectral features are attributed to Fabry-Perot-type fringes. The pump fre- 
quencies (shaded box) were tuned at the red edge through half the electronic 
band gap E G = 1.12 eV of silicon (vertical dashed line), and the probe frequen- 
cies at the blue edge of the stopband. The dashed curve is a calculation with 
the Scalar Wave Approximation in the region of interest. ReflExtAllD4.eps 
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Fig. 3. Differential reffectivity AR/R versus probe delay At taken at different 
pump frequencies E Pump and at probe frequency -Ep ro bc = 1.13 eV, sample do- 
main Al. At At = ps, the pump and probe are coincident, and the differential 
reflectivity (Ai?/i?) coinc decreases. At At = 1 ps, the differential reflectivity 
(AR/ R)pc has increased due to the dispersion of the free carriers (FC). The 
peak pump intensity varies between Jp ump = 10 ± 1 GWcm _2 (_Ep ump = 0.516 
eV) and ip ump = 25 ± 2 GWcm _2 (_Ep ump = 0.75 eV), and the probe intensity 
was /probe = 3 ± 2 GWcm -2 . DiffdRRvst.eps 
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Fig. 4. (a): Measured linear reflectivity versus probe frequency (open squares) 
on domain D4 compared to the scalar wave approximation (solid curve), (b) 
The differential reflectivity (AR/R) FC at delay At = 1 ps caused by the free 
carriers shows mostly dispersive features, as seen from the symmetric variation 
of AR/R around 0. When increasing the pump frequency from 0.54 (solid 
triangles) to 0.62 eV (open triangles), these dispersive features increase in 
magnitude, (c): At At = ps, dispersive as well as absorptive features are 
observed in (Ai?/i?) coinc , increasingly so when increasing the pump frequency 
from E Pump = 0.54 eV to 0.62 eV. For E Pwap = 0.54 eV, we find An' = 
0.44 x 10" 3 and n" = 0.46 x 10~ 3 , while for E Pump = 0.62 eV the fits give 
An' = 1.1 x 10~ 3 and n" = 1.1 x 10" 2 . The pump intensity was Jp ump = 36 ±4 
GWcm- 2 (£ Pump = 0.54 eV) and 7 Pump == 46 ± 5 GWcm- 2 (£ Pump = 0.62 eV). 
AlltOpsAlltlps.eps 
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Fig. 5. Instantaneous changes of the complex refractive index, (a): Imaginary 
refractive index (n", solid squares) vs. pump frequency Ep ump (lower axis) ob- 
tained from fits in fig. HI The dashed curve is a fit of n" to an exponential. The 
right scale shows the corresponding absorption coefficient. For comparison, 
we plot the absorption coefficient for low pressure chemical vapor deposited 
p-Si (solid curve) annealed at 545°C, similar to the backbone of our wood- 
pile crystals. Upper abscissa is the sum of the pump and probe frequencies 
-Erotai- (b): Change in real refractive index An' (left scale, solid squares) vs. 
pump frequency -Ep ump from our nondegenerate measurements. The theoreti- 
cal relation predicted by [39] is also shown (solid curve, right scale). We plot 
degenerate measurements for comparison: Open squares from [13] and open 
circles from [15]. Data from [14] exceeds the scale by one order of magnitude. 
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Fig. 6. Nonlinear coefficient (3u vs. pump frequency (solid squares) obtained 
from data as in fig. HI Dashed curve is a calculation from Ref. [39] for c-Si. 
Other data are from [13] (c-Si, open squares), [15] (c-Si, open circles), [14] 
(c-Si, triangles) and [41] (p-Si, diamonds), beta.eps 
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Fig. 7. Differential reffectivity at coincidence vs. pump frequency extracted 
from figure [3] (domain Al) and corrected for pump beam parameters (squares, 
left scale). Circles are complex values of n extracted from fig. \5\ (domain 
D4) and reinserted in the extended SWA at domain Al. We have calculated 
the differential reflectivity expected from a nonlinear, purely dispersive direct 
bandgap material (solid curve) [39] by inserting the theoretical relation for 
ri2 into the extended SWA. The 'purely dispersive assumption' deviates from 
the data for pump frequencies above -Ep ump = 0.69 eV, or .&rotai — 1-82 eV. 
KerrScaledwTheory.eps 
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Fig. 8. Nondegenerate Figure of Merit (NFOM) for instantaneous switching 
vs. pump and probe frequency for polysilicon (see eq. ED- For E Purap < \E G 
(=0.56 eV), pump absorption is absent as we do not consider more than sec- 
ond order photon absorption. If in addition -Ep ump + -Ep ro bc < -^Opt = 1-5 
eV, the total absorption is leading to large NFOM outside of the scale of 
graph (crossed area). The black curve indicates the parameter space for which 
a maximum NFOM is reached if three photon absorption is neglected. At 
-Epump = 0.8 eV, n2(-Epum P ) crosses cm 2 GW _1 , and so the NFOM vanishes. 
For -Erotai >> -Eopt, the NFOM is < 1 x 10~ 5 (hatched area). The parameter 
space in our present measurements are bounded by the box (squares). The 
only other measurement of ri2 on a Si based photonic structure is [16] (circle). 
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Fig. 9. (a) Relative change in frequency AE/Eq for a given change in the rela- 
tive average refractive index An av /n av . For comparison, the change according 
to Bragg's law AE/Eq = — An av /n av . (b) Change in reflectivity for a given 
change in imaginary part of the refractive index n" . CmpEMMSWA.eps 
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Fig. 10. Reflectivity calculated with the scalar wave approximation, with and 
without the thin layer of SiN. CmpSiNnoSiN.eps 
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